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ORIGINAL  ARTICLE 

Molecular  indicators  of  stress-induced  neuroinflammation 
in  a  mouse  model  simulating  features  of  post-traumatic 
stress  disorder 

S  Muhie1'2,  A  Gautam3,  N  Chakraborty1,  A  Hoke1,  J  Meyer h off1,  R  Hammamieh3  and  M  Jett3 


A  social-stress  mouse  model  was  used  to  simulate  features  of  post-traumatic  stress  disorder  (PTSD).  The  model  involved  exposure  of 
an  intruder  (male  C57BL/6)  mouse  to  a  resident  aggressor  (male  SJL)  mouse  for  5  or  10  consecutive  days.  Transcriptome  changes  in 
brain  regions  (hippocampus,  amygdala,  medial  prefrontal  cortex  and  hemibrain),  blood  and  spleen  as  well  as  epigenome  changes 
in  the  hemibrain  were  assayed  after  1-  and  10-day  intervals  following  the  5-day  trauma  or  after  1-  and  42-day  intervals  following 
the  10-day  trauma.  Analyses  of  differentially  expressed  genes  (common  among  brain,  blood  and  spleen)  and  differentially 
methylated  promoter  regions  revealed  that  neurogenesis  and  synaptic  plasticity  pathways  were  activated  during  the  early 
responses  but  were  inhibited  after  the  later  post-trauma  intervals.  However,  inflammatory  pathways  were  activated  throughout  the 
observation  periods,  except  in  the  amygdala  in  which  they  were  inhibited  only  at  the  later  post-trauma  intervals.  Phenotypically, 
inhibition  of  neurogenesis  was  corroborated  by  impaired  Y-maze  behavioral  responses.  Sustained  neuroinflammation  appears  to 
drive  the  development  and  maintenance  of  behavioral  manifestations  of  PTSD,  potentially  via  its  inhibitory  effect  on  neurogenesis 
and  synaptic  plasticity.  By  contrast,  peripheral  inflammation  seems  to  be  directly  responsible  for  tissue  damage  underpinning 
somatic  comorbid  pathologies.  Identification  of  overlapping,  differentially  regulated  genes  and  pathways  between  blood  and  brain 
suggests  that  blood  could  be  a  useful  and  accessible  brain  surrogate  specimen  for  clinical  translation. 
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INTRODUCTION 

Stress-induced  maladaptive  responses  are  major  health  problems 
affecting  significant  proportions  of  emergency  responders, 
soldiers,  victims  of  abuse,  and  survivors  of  accidents,  violence 
and  natural  disasters.  Varied  types  and  durations  of  traumatic 
exposure  (single  or  episodic,  acute  or  chronic  stressful  events) 
have  been  shown  to  induce  post-traumatic  stress  disorder  (PTSD) 
in  human1,2  and  PTSD-like  phenotypes  in  animal  models.3-6 

Stress-induced  PTSD  is  a  complex,  poorly  understood  and 
growing  problem.  In  its  chronic  form,  PTSD  includes  multisystem 
disorders  with  comorbidities  encompassing  both  physical  and 
psychiatric  maladies.  Despite  being  prevalent,  robust  diagnoses 
and  prognoses  of  PTSD  are  limited  and  optimal  therapeutic 
targets  are  lacking.  The  primary  challenge  is  the  difficulty  of 
detecting  and  treating  early  stages  of  the  disorder,  as  well  as 
predicting  clinical  courses  and  outcomes.  Even  known  risk  factors 
such  as  genetic,  epigenetic,  endocrine  and  demographic  factors 
(sex,  education,  income  and  social  status),  prior  adverse  life 
events,  or  a  history  of  early  trauma,  do  not  effectively  predict 
disease  severity  or  chronicity.7  Moreover,  the  widely  practiced  self- 
reported  and  physician  opinion-based  diagnoses  lack  objectivity. 

Because  targets  for  early  intervention  remain  elusive,  a  better 
understanding  of  the  molecular  underpinnings  of  PTSD  patho¬ 
genesis  is  needed  to  meet  the  challenges  of  disease  prognosis, 
diagnosis  and  treatment.  This  includes  identifying  molecular 
alterations  in  signaling  pathways  and  processes  that  underlie 


the  response  to  PTSD  triggering  stress.  Differentially  altered 
pathways  and  processes,  in  conjunction  with  pre-existing  genetic 
and  epigenetic  factors,  are  likely  to  mediate  stress-induced 
pathologies,  and  contribute  to  the  onset,  development,  severity 
and  persistence  of  the  disorder. 

Increasing  evidence  highlights  the  involvement  of  persistent 
dysfunctions  of  immune  response  pathways  in  psychiatric  and 
psychopathological  diseases,8  including  PTSD.9  Particularly,  stress- 
induced  systemic  inflammation  has  been  associated  with  learning 
and  memory  impairments,  leading  to  cognitive  and  behavioral 
deficits.10  For  example,  stress-activated  nuclear  factor-xB  (NF-xB; 
and  interleukin  (IL)-1  (3)  signaling  has  been  reported  to  inhibit 
hippocampal  neurogenesis  potentially  leading  to  behavioral 
disorder.11  In  addition  to  transcription-mediated  modulation, 
long-lasting  molecular  mechanisms  such  as  epigenetic  modifica¬ 
tions  have  been  implicated  in  mediating  the  persistent  impact  of 
stress-induced  inflammation  on  cognition  and  behavior.10  Despite 
these  promising  findings,  most  studies  attempting  to  connect 
the  effect  of  inflammatory  molecules  with  cognition  and 
behavioral  deficits  are  limited  by  their  focus  on  single  time  point 
observations. 

A  temporal  delineation  of  PTSD  development  from  the  nuclear 
stressor  event  to  altered  molecular  events  is  critical  for  identifying 
prognostic  markers.  However,  availabilities  of  molecular  data 
directly  connecting  clinical  phenotypes  to  signaling  pathways  and 
cellular  processes  are  currently  limited.  Furthermore,  molecular 
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mechanisms  and  outcome  predictors  determining  the  temporal 
trajectories  of  stress-induced  pathologies  leading  to  PTSD  are  yet 
to  be  elucidated.  In  this  study,  we  examined  longitudinal  multi- 
omics  alterations  across  tissues  and  behaviorally  relevant  brain 
regions  and  blood  from  a  social-stress  mouse  model  simulating 
features  of  PTSD.  Such  approaches  have  the  potential  to  provide  a 
better  understanding  of  molecular  alterations  leading  to  PTSD. 
First,  temporal  trajectories  should  help  to  delineate  the  type  and 
timing  of  molecular  events  underlying  adaptive  versus  maladap¬ 
tive  reactions  to  traumatic  stress.  Second,  a  shift  in  focus  from  a 
single  cytokine  or  neuronal  factor  to  networks  of  molecular 
interactions  and  their  predicted  activation  status  should  provide 
molecular  mechanisms  contributing  to  PTSD,  as  well  as  the 
molecular  milieu  in  which  comorbidities  are  maintained. 

The  social-stress  model,  used  in  this  study,  involved  random 
exposure  of  an  intruder  male  C57BL/6  mouse  to  a  trained 
aggressor  resident  male  SJL  mouse.  The  subject  C57BL/6  mouse 
was  later  exposed  to  trauma  reminders  (aggressor  barrier 
avoidance  tests)  at  various  post-trauma  intervals.12  Such  random 
exposure  to  a  conspecific  aggressor  was  found  to  elicit  important 
behavioral  features  reported  to  occur  in  PTSD  patients  such  as 
avoidance  of  trauma  reminders  (marked  avoidance  of  aggressor 
barrier),  jumping  during  handling  (increased  vigilance)  or  freezing 
(social  withdrawal)  and  impaired  cognition.12-14  Using  this  model, 
we  profiled  transcriptome  alterations  in  brain  regions  (hippocam¬ 
pus,  amygdala,  medial  prefrontal  cortex  and  hemibrain),  blood 
and  spleen  as  well  as  DNA  methylome  changes  in  the  hemibrain 
at  multiple  time  points. 

Transcripts  with  overlapping  expression  profiles  among  brain, 
blood  and  spleen  showed  activated  inflammation,  as  well  as 
inhibited  neurogenesis  and  impaired  synaptic  plasticity  at  longer 
post-trauma  delays.  Activation  status  of  signaling  pathways  at 
multiple  time  points  demonstrated  the  systemic  effect  of 
inflammation  on  neurogenesis  and  on  cognition  (learning  and 
memory),  while  also  accounting  for  the  chronic  pain  and  physical 
complaints  of  PTSD  patients.  Identification  of  overlapping 
differentially  expressed  genes  (DEGs)  and  pathways  between 
blood  and  brain  suggests  that  blood  can  be  used  as  an  accessible 
brain  surrogate  sample  for  clinical  translation. 

Our  findings  suggest  that  stress-activated  inflammatory 
responses  play  a  critical  role  in  the  development  and  maintenance 
of  PTSD-like  disorders  and  comorbid  conditions.  Hence,  diagnosis 
and  treatment  to  alleviate  stress-induced  PTSD  might  eventually 
include  strategies  for  addressing  the  inflammatory  problem  to 
significantly  improve  patients’  quality  of  life. 

MATERIALS  AND  METHODS 

Mouse  model  and  behavioral  evaluations 

Details  about  the  mouse  model  and  PTSD-related  behavioral  characteriza¬ 
tions  of  the  model  have  been  reported  earlier.12,13  Briefly,  naive  C57BL/6 
male  mice  (8-10  weeks  old)  were  single-housed,  and  mice  (A/ =  5  per 
group)  were  randomly  assigned  to  aggressor-exposed  (Agg-E)  and  cage- 
control  (C-ctrl)  groups  for  each  time  point.  Each  Agg-E  mouse  was  kept  in  a 
mesh  box  inside  the  cage  of  a  trained,  highly  aggressive,  resident  SJL 
mouse  for  6  h/day  for  either  5  (T5)  or  10  (T10)  consecutive  days.  Each  day, 
a  different  aggressor  mouse  was  used,  with  the  Agg-E  mouse  subjected  to 
two  to  three  direct  exposures  to  the  aggressor  mouse,  at  random,  for  1  min 
or  1 0  strikes,  whichever  came  first.  This  set-up  added  both  uncontrollability 
and  randomness  factors  to  the  social  defeat  stress  (to  simulate  some 
aspects  of  PTSD).  The  defeated  mice  were  inspected  before  and  after  each 
defeat  session  for  scratches  or  wounding  to  exclude  any  mouse  with  skin- 
deep  wounds  from  the  experiment.  Control  (never  exposed  to  aggressor) 
mice  were  identically  housed  for  the  6-h  period  in  a  separate  room  (for  the 
same  number  of  days).  An  hour  before  cervical  dislocation,  each  mouse 
was  tested  for  avoidance  behavior  by  measuring  the  time  spent  in  an 
‘interaction  zone’  around  a  caged  aggressor,  physically  separated  by  a 
plastic  divider  with  holes,  which  allowed  visual,  olfactory  and  auditory 
contact  during  the  avoidance  behavioral  assessment.  Summaries  of 


experimental  conditions  and  collected  tissues  are  shown  in  Supplemen¬ 
tary  Figure  1  and  Supplementary  Table  1.  Experiments  and  procedures  on 
mice  were  approved  by  the  local  (Walter  Reed  Army  Institute  of  Research) 
Institutional  Animal  Care  and  Use  Committees,  and  were  conducted  in 
strict  accordance  with  AAALAC  Guidelines. 

Y-maze  test 

The  spontaneous  alternation  task  was  performed  using  the  Y-maze 
apparatus  with  three  identical  arms  that  were  12.7  cm  wide,  20.32  cm  long 
and  evenly  separated  at  120°  from  each  other  (which  were  arbitrary  called 
arms  A,  B  and  C).  The  Y-maze  test  was  performed  1  day  after  T1 0  Agg-E  (or 
control),  and  behaviors  were  observed  for  7  min,  following  1-min  acclima¬ 
tion  to  the  maze.  Combinations  of  three  unique  arms  such  as  ‘ABC’,  ‘BAC’ 
and  ‘BCA’  were  used  as  spatial  parameters,  and  percent-alternation  rates 
were  calculated  using  the  number  of  possible  alternations: 

„ .  r  ,  .  number  of  alternations  (combination  of  unique  arms)  x  100 

%  of  alternation  rate  = , - , - 7 — ; - 7 - . - . - - - - r 

total  number  of  possible  alternations  (total  number  of  possible  arms) 


Sample  collection 

Mice  were  killed  by  cervical  dislocation,  alternating  between  Agg-E  and 
C-ctrl  mice  to  control  for  time  of  day  effects.  Blood  from  each  mouse  was 
collected  by  percutaneous  cardiac  puncture  and  stabilized  with  the 
PAXgene  Blood  RNA  solution  (Qiagen,  Valencia,  CA,  USA).  The  spleens  were 
also  collected  to  identify  more  stable  (less  transient)  overlapping  molecular 
alterations  with  (the  more  dynamic)  blood.  Brain  regions,  hemibrain, 
hippocampus  (HC),  amygdala  (AY)  and  medial  prefrontal  cortex  (MPFC), 
were  dissected  following  a  previously  reported  procedure.14 

RNA  isolation 

Total  RNA  was  isolated  from  homogenized  brain  regions  and  spleen  in  the 
TRIzol  reagent  (Invitrogen,  Thermo  Fischer  Scientific,  Carlsbad,  CA,  USA) 
and  from  blood  using  the  PreAnalytiX  PAXgene  Blood  RNA  kit  (Qiagen) 
following  the  manufacturers’  protocols.  RNA  was  evaluated  using  the 
NanoDrop  spectrophotometer  (Thermo  Fischer  Scientific)  and  Tapesta- 
tion2200  (Agilent  Technologies,  Santa  Clara,  CA,  USA). 

Assays 

Expression  arrays,  DNA  methylation  (tiling)  arrays  and  validation  assays. 
Expression  microarrays  for  all  tissues  were  performed  using  Agilent’s 
Genome-wide  Mouse  Expression  arrays  (GE_4x44Kv2  two-color),  whereas 
DNA  methylation  tiling  arrays  for  hemibrain  were  carried  out  using 
Agilent’s  Mouse  CpG-lsland  arrays  (2x105K),  both  following  the 
manufacturer’s  protocol.  Hybridized  microarray  slides  were  scanned  using 
Agilent’s  G2505C-US09493743-Scanner,  and  images  were  processed  using 
the  default  set-up  of  Agilent’s  Feature  Extraction  Software  v10.7. 

Processed  data  sets  were  filtered  to  exclude  probes  with  missing  values 
in  more  than  one  sample,  and  quantile  normalized  using  Limma, 
R-Package  (www.bioconductor.org).  Normalized  data  sets  were  analyzed/ 
compared  for  stress  effects  using  the  moderated  T-test  of  the  Limma 
package. 

DEGs  from  blood  were  validated  by  measuring  plasma-protein  levels. 
Blood  samples  were  collected  in  50  pi  of  buffered  sodium  citrate  (0.105  m 
— 3.2%)  and  centrifuged  to  separate  plasma,  which  was  snap-frozen.  A 
measure  of  70  pi  of  plasma  was  shipped  to  Rules-Based-Medicine,  Austin, 
TX,  USA,  for  Rodent-MAPv.20  antigens  assay  using  the  Luminex 
technology. 

DEGs  from  hemibrain  were  validated  using  QuantiGene2.0  Multiplex- 
Bead-Assay  (QuantiGene,  Affymetrix,  Sanata  Clara,  CA,  USA),  whereas  DEGs 
from  spleen  were  validated  using  nanoString  transcript  assay  (nCounter 
GeneExpressionCode-Sets,  NanoString  Technologies,  Seattle,  WA,  USA). 

Inflammatory  panel  assays.  Mouse-Inflammation-Panel  v2  was  used  for 
AY,  HC  and  MPFC  samples  according  to  the  manufacturer’s  protocol 
(NanoString  Technologies).  Briefly,  total  RNA  (100  ng)  was  hybridized 
overnight,  followed  by  post-hybridization  processing.  Data  sets  were 
analyzed  using  the  recommended  settings  of  the  nSolver-pipeline 
(NanoString  Technologies)  for  background  subtraction,  quality  control  on 
samples/lanes  and  the  downstream  data  analyses  as  follows:  normal¬ 
ization,  statistical  significance  testing  and  fold  change  (FC)  calculations. 
Transcripts  that  passed  P<  0.1  and  FC>1.5  in  each  region  (AY,  HC  or 
MPFC)  at  one  or  more  time  points  (T5R1,  T10R1,  T5R10  and  T10R42  except 
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Figure  1.  (a)  Differentially  expressed  genes  (DEGs)  in  blood,  hemibrain  and  spleen  of  trauma-exposed  C57BL/6  male  mice  across  different  time 

points.  DEGs  (blue  and  red  dots)  from  each  tissue  at  each  time  point  were  identified  using  moderated  7-test  at  P  <  0.05  and  fold  change  (FC) 
>  1 .5  filters.  The  scatter  plots  summarized  the  number  and  directions  of  DEGs  from  each  tissue  at  each  time  point  including  their  significance 
levels  and  effect  sizes  (in  log  FCs).  The  gene  symbols,  P-values  and  log2  FCs  for  the  top  DEG  from  each  tissue  and  time  point  are  given  in 
Supplementary  Table  2.  (b)  Pathways  and  processes  related  to  inflammatory  responses  and  innate  immunity.  These  were  identified  using 
initially  unbiased  enrichments  for  overlapping  DEGs  across  tissues.  Predicted  activation  states  (purple  and  yellow  bars)  were  colored  based  on 
predicted  activation  or  inhibition  z-scores  of  pathways  or  processes  enriched  at  day  42  after  the  10-day  aggressor  exposure  session  (T10R42  ), 
which  possibly  corresponds  to  persistent  post-traumatic  stress  disorder  in  human.  The  identity,  FC  and  significance  levels  of  genes  associated 
with  each  of  these  pathways  and  process  are  given  in  Supplementary  Table  3.  (c)  Differentially  methylated  promoter  probes  in  hemibrain. 
Differentially  methylated  promoter  CpG  islands  (blue  and  red  dots)  were  identified  using  moderated  7-test,  P  <  0.05,  FC>  1.5.  The  scatter  plots 
showed  promoter  regions  of  more  genes  were  differentially  methylated  due  to  the  longer  trauma  session  (10-day  aggressor  exposure 
session),  and  even  more  so  at  the  longer  post  aggressor  exposure  days  (at  T10R42).  (d)  Processes  and  pathways  associated  with  differentially 
methylated  probes  from  hemibrain  at  the  promoter  region  {q  <  0.05).  Predicted  activation  or  inhibition  states  of  these  pathways  and 
processes  were  colored  based  on  the  activated  or  inhibited  methylation  z-scores  of  promoter  CpGs  of  genes  significantly  associated  with  the 
corresponding  pathways  and  processes.  Promoter  regions  of  genes  significantly  associated  with  inflammatory  pathways  and  processes  were 
predicted  to  have  inhibited  DNA  methylation  states,  whereas  promoter  regions  of  genes  significantly  associated  with  neurogenesis  were 
predicted  to  have  activated  DNA  methylation  pattern.  (These  pathways  and  processes  were  enriched  using  hypergeometric  test  with  false 
discovery  rate  of  q  <  0.05.) 


3 


MPFC-T10R42,  due  to  sample  availability)  were  used  for  further  analyses 
and  downstream  visualization. 

Data  analysis 

Differentially  methylated  promoter  regions  and  DEGs,  as  determined 
between  Agg-E  and  C-ctrl  mice,  were  identified  at  the  four  time  points 
using  the  Moderated  T-test  in  Limma  (false  discovery  rate,  q^  0.05)  and 
0.05,  respectively).  Time  effects  and  trends  in  transcriptome  changes 
across  post-trauma  days  (for  overlapping  DEGs  among  brain,  blood  and 
spleen)  were  identified  using  Bayesian  and  temporal  clustering  methods.15 

Hypergeometric  test  {q  <  0.05),  Bingov3.0.3  (http:  apps.cytoscape.org/ 
apps/bingo)  and  ClueGOv2.2.5  (http://www.cytoscape.org),  Fisher’s  exact 
test  (Ingenuity,  Redwood,  CA,  USA),  Network  Analyst16  and  Gene  Set 
Enrichment  Analysis  (http://www.broadinstitute.org/gsea)  were  used  to 
determine  an  unbiased  gene  ontology  and  pathway  enrichments  of 
overlapping  DEGs.  Significant  processes  and  pathways  were  graphed 
and/or  visualized  using  R-programmingv3.3.0,  Cytoscapev3.4  (https://cran. 
r-project.org/),  GephivO.9.1  (https://gephi.org)  and  Ingenuity  Pathways 
Analysis. 


Activation  z-scores  for  significantly  enriched  pathways  and  processes 
were  calculated  to  predict  activation  states  based  on  intermolecular 
connectors  (edges).17,18  Experimental  directions  of  gene  expression  for 
members  of  the  pathway  or  the  network  were  compared  with  literature- 
derived  activation  directions  (as  curated  in  Ingenuity  Pathways  Analysis) 
for  either  activating  or  inhibiting  effects  (cause-effect  relationship  and 
direction  of  effect).17,18 

Functional  networks  and  bar  graphs  for  significant  pathways  and 
processes  indicating  the  directions  and  number  of  up-  or  downregulated 
DEGs,  and  activation  state  of  the  corresponding  pathway  or  process 
(values  and  directions  of  z-scores),  were  constructed  using  Cytoscapev3.4.0 
and  ggplot2  package  of  R,  and  visualized  using  GephivO.9.1  and 
RstudiovO.99.484  (https://www.rstudio.com). 

RESULTS 

DEGs  and  methylated  promoter  regions 

DEGs  from  each  tissue  at  each  time  were  identified  using  a 
moderated  f-test  (at  P<  0.05,  and  FC>1.5).  DEGs  that  passed 
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Figure  2.  (a)  Functional  network  for  differentially  expressed  genes  common  among  blood,  spleen  and  hemibrain  that  were  associated  with 

neurological  processes  and  pathways.  Gene  node  were  colored  based  on  the  expression  values  at  day  42  after  the  10-day  aggressor  exposure 
session  (T10R42),  and  pathway  nodes  were  colored  based  on  their  predicted  activation  states  at  later  time  points  (T5R10  and  T10R42).  (b) 
Functional  network  for  processes  and  pathways  associated  with  differentially  methylated  probes  from  hemibrain  at  the  promoter  region. 
Gene  nodes  were  colored  based  on  the  methylation  values  of  promoter  CpG  islands  of  the  gene  at  day  42  after  the  1 0-day  aggressor  exposure 
session  (T10R42),  and  pathway  nodes  were  colored  based  on  predicted  activation  of  methylation  states  of  the  pathways  as  inferred  from 
methylation  status  of  promoter  regions  of  genes  involved  in  the  respective  pathway.  Hence,  in  general,  the  actual  activation  of  the  pathway  is 
expected  to  be  opposite  of  what  is  shown  here.  (These  pathways  and  processes  were  enriched  using  hypergeometric  test  with  false  discovery 
rate  of  q  <  0.05.) 


P<  0.05  and  FC>1.5  filter,  and  were  overlapping  among  blood, 
hemibrain  and  spleen  (Figure  la),  showed  two  different  temporal 
trends:  1874  DEGs  increased  in  expression,  whereas  1757  DEGs 
decreased  in  expression  with  increasing  post-trauma  periods 
(Supplementary  Figure  2).  These  sets  of  DEGs  were  mainly 
associated  with  immune  response  (innate  immunity  and  inflam¬ 
matory  responses;  Figure  1b)  and  neurological  system  processes 
(mainly  neurogenesis;  Figure  2a). 

At  later  post-trauma  days  (T5R10  and  T10R42),  pathways  related 
to  innate  immunity  and  inflammatory  responses,  such  as 
interferons,  Toll-like  receptors  and  tumor  necrosis  factor-a  (TNF- 
a),  and  interleukin  signaling,  chemotaxis  and  NF-KB-regulated 
transcriptions  were  activated.  On  the  contrary,  signaling  pathways 
implicated  in  anti-inflammatory  responses  (SOCS1,  SOCS3,  TOB1, 
CTLA4,  IL-1RN  and  IL-10  signaling)  were  inhibited  (Figure  1b; 
Supplementary  Table  3).  Interestingly,  however,  pathways  related 
to  CD3,  CD4,  CTLA4,  TOB1,  IDOI  and  TGFB  signaling  were 
suppressed,  suggesting  inhibited  T-helper  cell-type-2  immunity 
(Figure  1b;  Supplementary  Table  3),  whereas  most  immune 
responses  were  activated.  Biological  processes  related  to  neuro¬ 
genesis,  synaptic  plasticity,  learning  or  memory,  mTOR,  NOTCH 
and  calcium  signaling  were  inhibited,  whereas  apoptosis  of 
neurons,  WNT-  and  cyclic  AMP  signaling  was  activated  at  later 
post-trauma  periods  (Figure  2a). 

Consistent  with  the  above  findings,  upregulated  central  nodes 
in  the  regulatory  networks  for  the  overlapping  DEGs  (Irf8,  Ddx58, 
NfKbl,  116,  Irfl,  Cebpb,  Bcl2  and  Trim8)  were  involved  in 
inflammatory  responses  (Supplementary  Figure  3),  whereas 
downregulated  central  nodes  were  involved  in  immunodeficiency 
processes  (Pik3r1)  and  modulating  synaptic  plasticity  (Tall).19 

Enrichments  of  significantly  (q  <  0.05)  hyper-  and 
hypo-methylated  promoter  regions  in  the  hemibrains  of  the 


longer-period  (T10)  stress  group  at  the  later  (T10R42)  post-trauma 
time  point  (Figure  1c)  showed  inhibited  promoter-methylation 
patterns  of  genes  related  to  inflammation  and  innate  immune 
responses,  and  increased  methylation  patterns  at  promoter 
regions  related  to  neurogenesis,  synaptic  plasticity,  activities  of 
dopaminergic  synapses,  circadian  rhythm,  telomere  maintenance, 
and  FOS  and  mTOR  signaling  (Figures  Id  and  2b).  This  is 
consistent  with  enrichments  and  activation  states  of  DEGs 
common  among  brain,  blood  and  spleen  (Figures  1b  and  2a). 

In  summary,  at  42  days  after  the  10  days  of  trauma  (T10R42) — 
believed  to  roughly  approximate  chronic  PTSD — functional  net¬ 
works  of  DEGs  common  among  brain,  blood  and  spleen  showed 
inhibition  of  neurogenesis,  synaptic  plasticity  and  processes 
related  to  learning  and  memory,  which  was  further  supported 
by  changes  in  DNA  methylation  patterns  in  the  hemibrain.  On  the 
other  hand,  signaling  pathways  related  to  inflammatory  responses 
were  activated,  whereas  anti-inflammatory  pathways  and  T-helper 
cell-type-2  immunity  were  inhibited. 

Additional  molecular  changes  that  were  associated  with  PTSD- 
comorbid  pathologies  (metabolic  disorders,  insulin  signaling 
disruption,  oxidative  stress,  dysregulation  of  mitochondrial  and 
telomere  maintenance)  were  also  significantly  associated  with 
commonly  altered  transcripts  (Table  1;  Supplementary  Table  4). 

Molecular  alterations  in  brain  regions  implicated  in  fear  and 
traumatic  stress  responses 

We  further  investigated  stress-induced  molecular  alterations  in 
brain  regions  implicated  in  traumatic  stress  response:  AY,  HC  and 
MPFC.  Others  have  shown  that,  under  traumatic  stress,  the 
dentate  gyrus  of  HC  (an  active-site  of  neurogenesis  in  the  adult 
brain)  is  highly  affected  by  hippocampal  inflammation,20  amyg¬ 
dala  activity21  and  MPFC  neuronal  inputs,22  making  these  three 
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Table  1.  Differentially  regulated  genes  (common  among  blood,  hemibrain  and  spleen)  that  were  significantly  associated  with  signaling  pathways 
implicated  in  PTSD  comorbidities 


Signaling  pathway 


Differentially  expressed  genes 


Calcium  signaling  pathway 


cAMP  signaling  pathway 


Insulin  signaling  pathway 


Mitochondrial  biogenesis  and 

translation 

mTOR  signaling 

Oxidative  stress  response 

Telomerase  maintenance  and  aging 
Type  I  and  II  diabetes  mellitus 
VEGF  signaling  pathway 


WNT  signaling  pathway 


RYR2,  GNA14,  CYSLTR1,  ADORA2A,  ADRA1B,  PTGER1,  HRH2,  CALML5,  ATP2B1,  ATP2B2,  CD38,  ATP2A2, 
ERBB3,GNAL,  GNAS,  PDE1C,  PDE1A,  EGFR,  CACNA1S,  ADCY1,  HTR6,  EDNRB,  PLCB2,  ADCY3,  ADCY5, 
CAMK2B,  GRB2 

RYR2,  RAC2,  EP300,  ADORA2A,  CALML5,  ATP2B1,  ATP2B2,  ATP2A2,  ADCYAP1R1,  PIK3R1,  GIPR,  GNAS, 
PDE4D,  VAV3,  VAV2,  CACNA1S,  ACOX1,  ADCY1,  ATP1B1,  HTR6,  NFKB1,  PPP1CC,  AKT1,  ADCY4,  ADCY5, 
ADCY9,  ADORA3,  ADRA2A,  ADRBK1,  AKAP6,  CAMK2B,  CHRM1,  FLNA,  GLP2R,  GNAL,  GRM6,  GRM7,  OPRD1, 
OPRM1,  PDE1A,  PTGER4,  RAPGEF4,  RGS4,  RIMS2,  RLF,  RPS6KA1,  S1PR3,  GPR44 
PIK3CD,  ACACA,  GSK3B,  EIF4EBP1,  CALML5,  PDPK1,  PRKAB2,  PRKAB1,  GRB2,  SHC2,  PIK3R1,  PTPRF, 
PPP1R3E,  RPS6,  SOS1,  TSC1,  CBLC,  PPP1CC,  PTPN1,  EIF4E2,  AKT1,  INPP5D,  PTPRA,  ADIPOR1,  CD38,  FFAR1, 
GIPR,  GJA1,  I  RSI,  JAK2,  LARS,  MGEA5,  PFKFB2,  SERP1,  SOCS1,  TCF4,  TCF7L2,  TRH,  ADCY1,  ADCYAP1R1, 
ATP1B1,  CDK2,  CPLX1,  CPLX3,  CRKL,  CSNK2A1,  DOK5,  DOK7,  FGF23,  FOXOI,  GNAS,  I  KB  KB,  IRS4,  KIF5B, 
KLB,  LPIN1,  MY05A,  NRAS,  PTPN2,  RAF1,  RAPGEF1,  RAPGEF4RIMS2,  RYR2,  SHC4,  SMARCC1,  SOCS4,  STAT1, 
STXBP4,  CALM1,  CALM3,  CALM4,  CRK,  FGFR4,  IRS3,  KRAS,  MAP3K14,  MAP3K8,  MAP4K3,  MAP4K5,  MAPK6, 
PDE3A,  PIK3C2G,  PRKCB,  PRKCH,  PRKCQ,  RHOJ,  SNAP23,  VAMP2 

MEF2D,  HDAC3,  TFAM,  PPARGC1 B,  TGS1,  PRKAB2,  PRKAB1,  HELZ2,  MAPK12,  BAK1,  DAP3,  GFM2,  MRPL23, 
MRPL24,  MRPL33,  MRPL44,  MRPL49,  MRPL9,  MRPS15,  MRPS17,  MRPS23,  MRPS25,  MRPS33 
AKT1,  CDK2,  CLIP1,  EIF3D,  EIF4A1,  EIF4B,  EIF4E2,  EIF4EBP1,  IGF1,  PDPK1,  PIK3CD,  PIK3R1,  PPP2R5E,  RHOF, 
RPS3,  RPS6,  RPS6KA1,  RRAGD,  TSC1,  VEGFA,  PRKCB 

DUSP10,DUSP3,DUSP7,BCL2,MYC,DDIT3,STAT1,MAPK12,CD24A,TPO,NCF2,CAT,CYCT,CYBB,DNAJB14, 

DNAJB4,DNAJB6,DNAJB9,DNAJC10,DNAJC11,DNAJC16,DNAJC17,DNAJC5,DNAJC9,SEC63 

TP53,  BCL2,  MYC,  POLR2A,  AKT1,  SMC6,  NSMCE2,  NFKB1,  PPP2R5E,  XRCC6,  CABIN1,  CDK2,  HIRA,  HMGA2 

CD80,  FAS,  IRF 1 ,  H2-AB1,  H2-EB1,  H2-K1,  H2-M10.2,  H2-T22,  TRAV3N-3,  AKT1,  PIK3CD,  PIK3R1 

RAC2,  ELMOI,  BAIAP2,  PTK2,  CYBB,  PDPK1,  SHC2,  PIK3R1,  NCF4,  CTNNA1,  DOCK1,  VAV3,  VAV2,  CDC42, 

SHB,  SH2D2A,  MAPK12,  AKT1,  CTNND1,  GRB2,  ITGA1,  ITGA4,  ITGB1,  KRAS,  PIK3C2A,  PLA2G4A,  PRKCB, 

SOS1,  VCAM1,  VEGFA,  RAF1,  PTPN2 

MYH7,  MYH6,  MYH8,  CDH16,  CSNK1D,  GSK3B,  GNA14,  EP300,  ARRB2,  HDAC3,  HDAC8,  PCDH1,  PCDH9, 
PCDH7,  PCDHB5,  PCDHB6,  PCDHB4,  CSNK1A1,  WNT6,  LEF1,  GNG5,  CTNNA1,  PCDH20,  TGFBR1,  PCDH10, 
FRAT1,  MYC,  GNB1,  PCDHB11,  PCDHB10,  PCDHB16,  CTNNAL1,  CDH9,  FAT3,  FZD9,  FZD3,  FZD6,  DKK1, 
WNT5A,  ACVR1C,  PLCB2,  TCF7L2,  WNT2B,  RAC2,  VRK2,  CDC42,  CSNK1E,  CSNK2A1,  CSNK2A2,  CSNK2B, 
CTNND1,  DDIT3,  FBXW2,  HHEX,  HOXB9,  LDB1,  MITF,  PIAS4,  PYG02,  RNF138,  RYR2SENP2,  SFRP4,  TCF3, 
TCF4,  TGFB1 II,  TLE2,  TLE3,  VANGL2,  WNT3,  WNT5B,  AKT1,  CAMK2B,  CUL1,  G2E3,  HECW2,  NFATC1,  NFATC3, 
PPP2R5B,  PPP2R5E,  PPP3CC,  PRICKLE1,  PRKCB,  RAC3,  ROCK1,  ROCK2,  WWP2,  DAB2,  DLG1,  FHL2,  HIPK2, 
RAF1,  RUNX2,  SUMOI,  TRP53 


Abbreviation:  PTSD,  post-traumatic  stress  disorder. 
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brain  regions  noteworthy23,24  for  additional  molecular  profiling. 
Differentially  regulated  transcripts  in  the  AY,  HC  and  MPFC  across 
four  time  points  were  significantly  associated  with  inflammatory 
pathways  (TNF-a,  IFN-y,  IL-6,  chemokine,  and  TLR  signaling  and 
NF-kB  transcriptional  activities),  as  well  as  other  important 
signaling  pathways  such  as  histone  modifications,  oxidative  stress, 
neuronal  growth  factors,  and  WNT,  NOTCH,  (3-catenin  and  cyclic 
AMP  signaling  pathways  (Figure  3a). 

In  the  AY,  signaling  pathways  of  inflammatory  responses, 
neurogenesis  and  synaptic  plasticity  were  activated  earlier  (a  day 
after  stress)  in  response  to  trauma,  and  inhibited  at  later  post¬ 
trauma  days  (42  days  after  stress)  (Figure  3a;  Supplementary 
Figure  4).  Suppressed  expression  of  DEGs  in  the  AY  at  later  post¬ 
trauma  days  (Figure  3a;  Supplementary  Figure  5a)  is  likely  related 
to  the  role  of  the  AY  as  an  early  fear-response  center 25  Notably,  as 
reduced  basolateral  amygdalar  activity  or  lesions  suppress  adult 
neurogenesis,  DEG  suppression  in  this  region  may  prevent 
selective  activation  of  newborn  neurons26  and  contribute  to 
decreased  hippocampal  neurogenesis. 

In  the  HC  and  MPFC,  inflammatory  signaling  pathways  were 
activated  at  earlier  time  points  (T5R1)  and  remain  activated  at  later 
post-trauma  days  (T10R42),  whereas  neurogenesis  and  synaptic 
plasticity  were  inhibited  at  T10R42  (Figure  4a).  Particularly,  DEGs 
from  HC  were  associated  with  signaling  pathways  of  VEGF,  IFN-y, 
IL-6,  TNF-a,  NF-kB,  Toll-like  receptors,  chemotaxis,  calcium,  notch 
and  circadian  rhythm  (Supplementary  Figure  5a).  In  addition  to 
evidence  of  activated  inflammatory  pathways,  expression  levels  of 
important  pro-inflammatory  molecules  (IFN-y,  IL-6,  TNF-a,  IL-1B, 
CRP,  CCR2  and  others)  were  increased  in  the  HC  and  MPFC  at 


earlier  and,  to  a  lesser  extent,  at  later  post-trauma  days  (Figure  3b; 
Supplementary  Figures  5b). 

Activated  inflammatory  responses,  inhibited  neurogenesis  and 
inhibited  synaptic  plasticity  in  HC  and  MPFC  are  consistent  with 
the  activation  of  inflammatory  pathways  and  inhibition  of 
neurogenesis  for  overlapping  DEGs  among  the  hemibrain,  blood 
and  spleen  observed  at  later  post-trauma  days  (Figure  4b).  They 
are  also  consistent  with  downregulated  neuronal  growth  factors 
(BDNF,  DLG4,  IGF1,  VEGFA,  ITGAL  and  EGFR)  in  overlapping  DEGs 
at  later  post-trauma  days  leading  to  inhibited  neurogenesis 
(Figure  4b). 

To  further  corroborate  our  molecular  indicators  of  inhibited 
neurogenesis  in  Agg-E  mice,  we  used  a  (Y-maze)  spatial  learning 
and  recognition  memory  (novelty  cognition)  test,  which  is  reported 
to  involve  newly  generated  neurons,27-29  although  the  association 
of  this  test  with  neurogenesis  is  inconsistent  across  studies.30  The 
stressed,  Agg-E  mice  showed  poorer  performance  than  the 
C-ctrl  mice  (P  <  0.09)  at  T10R1  (Figure  4c),  suggesting  impaired 
neurogenesis  in  the  dentate  gyrus  of  the  HC,  which  corroborates 
the  molecular  findings  regarding  inhibited  neurogenesis. 

In  summary,  we  observed  activated  inflammatory  signaling 
pathways  early  on  (T5R1)  in  the  HC  and  MPFC,  and  activation 
persisted  even  at  42  days  after  the  10  days  of  trauma  (T10R42)  in 
HC  and  MPFC.  In  both  the  HC  and  MPFC,  we  observed  not  only  the 
activation  of  downstream  signaling  pathways  of  pro-inflammatory 
molecules  but  also  increased  expression  levels  of  the  pro- 
inflammatory  molecules  themselves  at  T5R1  (including  those 
shown  to  inhibit  neurogenesis:  TNF-a,  CRP,  IL-6,  IFN-y,  IL-1B  and 
CCR2)  as  verified  using  an  independent  nanoString  assay 
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Figure  3.  (a)  Activation  states  of  inflammatory  response,  WNT  signaling,  insulin  signaling,  telomere  dysregulation  and  growth  factor-related 

signaling  pathways  in  amygdala  (AY),  hippocampus  (HC)  and  medial  prefrontal  cortex  (MPFC)  at  four  different  time  points  of  day  1  of  after  the 
5  days  of  trauma  (T5R1 )  and  1 0  days  of  trauma  (T1 0R1 ),  day  1 0  of  after  the  5  days  of  trauma  (T5R1 0),  and  day  42  of  after  the  1 0  days  of  trauma 
(T10R42).  (b)  Select  pro-inflammatory  molecules  from  ~  120  transcripts  assayed  using  nanoString's  255  pro-inflammatory  probes.  About  120 
transcripts  (Supplementary  Figure  5a)  passed  fold  change >1.5  and  a  significance  of  P<0.1  in  AY,  HC  and  MPFC  at  least  at  onetime  point  in 
each  brain  region. 
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Figure  4.  (a)  Functional  networks  of  differentially  regulated  genes  in  hippocampus  (top  panel)  and  medial  prefrontal  cortex  (lower  panel)  at 

day  1  of  after  the  5  days  of  trauma  (T5R1 )  and  1 0  days  of  trauma  (T1 0R1 ),  day  1 0  of  after  the  5  days  of  trauma  (T5R1 0),  and  day  42  of  after  the 
10  days  of  trauma  (T10R42).  (b)  Regulatory  networks  of  differentially  expressed  genes  common  among  blood,  hemibrain  and  spleen  at  later 
post-trauma  days  (T5R10  and  T10R42).  (c)  Y-maze  assessment  of  aggressor-exposed  (stressed)  and  control  C57BL6  mice  at  T10R1. 
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(Figure  4b).  Even  at  42  days  after  the  10  days  of  traumatic  stress, 
inflammation  remained  activated  in  the  HC  (though  to  a  lesser 
extent).  By  contrast  to  inflammatory  responses,  pathways  involved 
in  neuronal  growth  factors  (neurogenesis)  were  inhibited. 
Activated  inflammation  and  inhibited  neurogenesis  at  T10R42  in 
brain  regions  implicated  in  fear  memory  and  extinction  (HC  and 
MPFC)  again  corroborate  our  findings  from  overlapping  DEGs 
among  hemibrain,  blood  and  spleen. 

Blood  transcripts  versus  plasma  proteins 

We  observed  concordance  between  levels  of  plasma  proteins  at 
two  time  points  (T5R1  and  T10R1)  and  expression  levels 
(directions)  of  corresponding  transcripts  in  the  whole  blood, 
despite  the  different  cohorts  of  mice  used  for  plasma  and  whole 
blood  (transcriptome)  assays.  Eight  genes  (TPO,  IL-1 B,  MPO,  MMP9, 
VEGFA,  CRP,  TIMP1  and  APCS)  were  consistently  upregulated  in 
both  acute  (T5)  and  chronic  (T10)  trauma  sessions  at  the  levels  of 
both  protein  and  transcript  expression  (Supplementary  Figure  6a). 
Within  each  exposure  session,  13  out  of  a  total  of  16  genes  at 
T5R1 ,  and  1 1  out  of  1 6  genes  at  T1 0R1  were  expressed  in  the  same 
direction  at  both  the  protein  and  transcript  levels  in  each 
exposure  session  (Supplementary  Figure  6a).  Stress-induced 
alterations  of  these  proteins,  along  with  dysregulation  of  immune 
response  pathways,  may  underlie  the  maladaptive  responses  that 
can  potentially  progress  to  specific  behavioral  disorders  and 
somatic  pathologies. 

Validation  of  transcripts  from  hemibrain  and  spleen 
We  also  validated  DEGs  from  the  hemibrain.  Although  some  of  the 
assayed  transcripts  were  not  among  significant  DEGs,  we 
observed  a  significant  positive  correlation  mostly  in  direction 
and  in  many  cases  in  magnitude  of  FCs  (Supplementary 
Figure  6b).  For  DEGs  from  the  spleen,  the  nanoString  assay 
showed  a  directional  correlation  with  microarray  findings 
(Supplementary  Figure  6c). 

Overall,  our  mouse  model  made  it  feasible  for  temporal 
sampling  of  important  brain  regions  and  organs  for  system-level 
analysis  of  molecular  events  that  are  otherwise  inaccessible  from 
patients.  This  allowed  us  to  identify  differentially  expressed 
transcripts,  methylated  promoter  regions  across  multiple  time 
points  and  associated  signaling  pathways  that  reveal  potential 
molecular  mechanisms  for  stress-induced  disorders.  Molecular 
changes  at  earlier  time  points  provided  molecular  signatures  for 
disease  onset  and  progression,  whereas  those  at  later  time  points 
are  indicators  of  disease  persistence. 

DISCUSSION 

We  have  identified  temporal  and  overlapping  DEGs  across 
multiple  tissues,  and  associated  gene  interaction  networks  in  a 
mouse  model  simulating  aspects  of  PTSD.  Such  approaches 
provided  comprehensive  assessment  of  the  molecular  under¬ 
pinnings  of  the  major  aspects  of  this  debilitating  disorder. 
Specially,  determining  the  activation  status  of  significantly  altered 
pathways  or  processes  has  made  it  possible  to  infer  behavioral 
responses  (phenotypic  outcomes)  of  stress-induced  molecular 
alterations. 

Our  integrative  and  consensus  network  analyses  of  the  various 
blood,  brain  and  spleen  data  sets  suggest  that  stress-induced 
inflammation  may  be  negatively  affecting  brain  regions  respon¬ 
sible  for  cognition  and  traumatic  fear  memory  extinction  leading 
to  behavioral,12,14  and  systemic  disorders  including  metabolic 
imbalance,31  further  dysregulation  of  immune  response,  disrup¬ 
tion  of  insulin  signaling,  dysfunction  of  mitochondrial  and 
telomeric  maintenances,  and  cardiovascular  problems.32  Our 
findings  suggest  that  activated  inflammation  at  earlier  time  points 
in  brain  regions  implicated  in  fear  memory  and  extinction  (HC  and 


MPFC)  might  cause  inhibited  neurogenesis  at  later  time  points. 
This  is  consistent  with  studies  that  show  causal  relationships 
between  hippocampal  inflammation  and  impaired  neurogenesis 
in  the  dentate  gyrus  20 

Normally,  chemokines,  cytokines  and  the  family  of  NF-kB  factors 
are  widely  expressed  in  zones  of  active  neurogenesis,  where  they 
are  critical  for  maintaining  brain  health  and  function.33,34  This 
includes  long-term  potentiation  of  hippocampal  CA1  region,35 
neurogenesis34,36  and  synaptic  plasticity  that  underlies  learning 
and  memory.37  However,  over-secretion  of  pro-inflammatory 
cytokines  has  been  associated  with  neurodegenerative  and 
neuropsychiatric  diseases.38,39  For  example,  elevated  levels  of 
IL-6,  TNF-a,  IFN-y  and  NF-kB  underlie  inflammation-induced 
cognitive  decline,40  C-reactive  protein-mediated  inflammation 
has  been  associated  with  PTSD  symptoms,40  and  IL-1  (3-induced 
chronic  neuroinflammation  impairs  adult  HC  neurogenesis41 
disrupting  production,  distribution  and  recruitment  of  new 
neurons  into  relevant  neural  networks  42  Even  mild  inflammation 
in  the  HC  is  linked  to  mood  changes,  memory  deficits,  reduced 
neurogenesis  and  partial  loss  of  principal  neurons,43  whereas 
chronic  peripheral  inflammation  has  been  shown  to  inhibit 
neurogenesis.44  Other  comorbidities  such  as  diabetes,  which  we 
infer  to  be  triggered  by  inflammation,  have  also  been  implicated 
in  impaired  hippocampal  neurogenesis.45  Moreover,  increased 
plasma  levels  of  pro-inflammatory  cytokines,  upregulated  targets 
of  NF-KB/Rel  transcription  factors46  and  molecular  indicators  of 
enhanced  inflammation  have  been  linked  to  PTSD  vulnerability  47 
and  to  comorbid  somatic  diseases  in  childhood-abuse48  and 
combat-related  PTSD.49 

Here  observation  of  prolonged  inflammation  and  inhibited 
methylation  of  promoter  regions  for  pro-inflammatory  molecules 
are  also  consistent  with  multiple  reports  that  show  a  correlation 
between  stress-induced  inflammation  and  chronic  pain  in 
PTSD.40,50-53  Stress-induced  chronic  inflammation  is  also  impli¬ 
cated  in  chronic  pain,53-55  which  may  be  responsible  for  the 
physical  complaints  possibly  contributing  to  the  impaired  social 
and  emotional  functions  of  PTSD  patients.56  Additional  long-term 
effects,  most  notably  epigenetic  changes  in  the  adult  brain,  could 
further  exacerbate  alterations  to  long-term  memory  and  synaptic 
plasticity,  contributing  to  PTSD-related  inhibition  of  neurogenesis 
and  hence  cognitive  deficit.  These  molecular  observations  were 
phenotypically  corroborated  by  the  impaired  Y-maze  performance 
of  the  stressed  mice,  a  behavioral  deficit  observed  for  we  believe 
the  first  time  in  an  animal  model  simulating  aspects  of  PTSD. 

These  transcriptome-  and  epigenome-based  correlative  findings 
(and  inferred  causal  relationships  between  inflammation  and 
PTSD-associated  pathologies)  support  the  notion  that  stressed- 
induced  inflammation  triggers  or  maintains  cascades  of 
comorbidities.57-59  In  short,  activated  neuroinflammation  is 
implicated  in  inhibited  neurogenesis  and  synaptic  plasticity 
(cognitive  deficits)  leading  to  the  development  and  persistence 
of  the  behavioral  manifestations  of  PTSD.  Whereas,  somatic 
inflammation  seems  to  be  directly  involved  in  tissue  damage, 
triggering  and  maintaining  somatic  comorbid  pathologies. 

The  association  between  PTSD  and  chronic  inflammation 
suggests  that  treatment  of  trauma-induced  inflammation  may 
have  wide-reaching  implications  for  ameliorating  PTSD  symptoms. 
Despite  inconsistencies  regarding  roles  of  specific  pro- 
inflammatory  molecules  in  PTSD,60  which  might  be  attributed  to 
psychotropic  medications,  variations  in  assay  methods  or  blood 
collection  time  of  the  day,46  inflammatory  pathways  may  serve  as 
therapeutic  targets.  Hence,  diagnosis  and  treatment  strategies 
designed  to  address  the  burden  of  PTSD  may  include  addressing 
the  inflammatory  problem  (immune  response  imbalances)  to 
alleviate  chronic  pain,  and  to  reverse  tissue  damage  and  cognitive 
deficits  potentially  rooted  in  inflammation. 

Identification  of  overlapping  molecular  changes  in  multiple 
tissues,61  particularly  the  presence  of  common  molecular 
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alterations  between  blood  and  brain,  might  confer  confidence  in 
using  blood  as  an  accessible  brain  surrogate  sample  for  clinical 
translation.  Such  overlapping  molecular  events  potentially  provide 
a  convenient  method  for  patient  monitoring  and  PTSD  evaluation, 
leading  to  better  prognosis,  diagnosis  and  treatment  strategies. 

Data  depositions 

Genome-wide  data  sets  reported  in  this  manuscript  are  deposited 
at  the  Gene  Expression  Omnibus  of  the  National  Center  for 
Biotechnology  Information. 

GSE85495:  http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? 

acc  =  GSE85495 

GSE68076:  http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? 

acc  =  GSE68076 

GSE45035:  http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? 

acc  =  GSE45035 

CONFLICT  OF  INTEREST 

The  authors  declare  no  conflict  of  interest. 

ACKNOWLEDGMENTS 

We  thank  Stacy-Ann  Miller,  Bintu  Sowe  and  Seshamalini  Srinivasan  for  nucleic  acid 
isolation  and  microarray  experiments;  Ross  Campbell,  Sharon  Bewick,  John  Clifford, 
Derese  Getnet  and  Duncan  Donohue  for  editing  the  manuscript;  and  US  Army 
Medical  Research  and  Materials  Command  (USAMRMC#:  09284002)  for  grant  support. 

DISCLAIMER 

Research  was  conducted  in  compliance  with  the  Animal  Welfare 
Act,  and  all  other  Federal  requirements.  The  views  expressed  are 
the  authors’  and  do  not  constitute  endorsement  by  the  US  Army 

REFERENCES 

1  O'Donnell  ML,  Varker  T,  Creamer  M,  Fletcher  S,  McFarlane  AC,  Silove  D  et  al. 
Exploration  of  delayed-onset  posttraumatic  stress  disorder  after  severe  injury. 
Psychosom  Med  2013;  75:  68-75. 

2  Van  Loey  NE,  van  de  Schoot  R,  Faber  AW.  Posttraumatic  stress  symptoms  after 
exposure  to  two  fire  disasters:  comparative  study.  PLoS  ONE  2012;  7:  e41532. 

3  Jovanovic  T,  Sakoman  AJ,  Kozaric-Kovacic  D,  Mestrovic  AH,  Duncan  EJ,  Davis  M 
et  al.  Acute  stress  disorder  versus  chronic  posttraumatic  stress  disorder:  inhibition 
of  fear  as  a  function  of  time  since  trauma.  Depress  Anxiety  2013;  30:  217-224. 

4  Reich  CG,  Iskander  AN,  Weiss  MS.  Cannabinoid  modulation  of  chronic  mild  stress- 
induced  selective  enhancement  of  trace  fear  conditioning  in  adolescent  rats. 
J  Psychopharmacol  2013;  27:  947-955. 

5  Daskalakis  NP,  Cohen  H,  Cai  G,  Buxbaum  JD,  Yehuda  R.  Expression  profiling 
associates  blood  and  brain  glucocorticoid  receptor  signaling  with  trauma-related 
individual  differences  in  both  sexes.  Proc  Natl  Acad  Sci  USA  2014;  111: 
13529-13534. 

6  Ronzoni  G,  Del  Arco  A,  Mora  F,  Segovia  G.  Enhanced  noradrenergic  activity  in  the 
amygdala  contributes  to  hyperarousal  in  an  animal  model  of  PTSD.  Psychoneur¬ 
oendocrinology  2016;  70:  1-9. 

7  Xue  C,  Ge  Y,  Tang  BH,  Liu  Y,  Kang  P,  Wang  M  et  al.  A  meta-analysis  of  risk  factors 
for  combat-related  PTSD  among  military  personnel  and  veterans.  PLoS  ONE  2015; 
10:  eOl  20270. 

8  Ma  L,  Piirainen  S,  Kulesskaya  N,  Rauvala  H,  Tian  L.  Association  of  brain  immune 
genes  with  social  behavior  of  inbred  mouse  strains.  J  Neuroinflammation  201 5;  1 2: 
75. 

9  Segman  RH,  Shefi  N,  Goltser-Dubner  T,  Friedman  N,  Kaminski  N,  Shalev  AY. 
Peripheral  blood  mononuclear  cell  gene  expression  profiles  identify  emergent 
post-traumatic  stress  disorder  among  trauma  survivors.  Mol  Psychiatry  2005;  10: 
500-513,  425. 

10  Donzis  EJ,  Tronson  NC.  Modulation  of  learning  and  memory  by  cytokines: 
Signaling  mechanisms  and  long  term  consequences.  Neurobiol  Learn  Mem  2014; 
115:  68-77. 

1 1  Koo  JW,  Russo  SJ,  Ferguson  D,  Nestler  EJ,  Duman  RS.  Nuclear  factor-kappaB  is  a 
critical  mediator  of  stress-impaired  neurogenesis  and  depressive  behavior.  Proc 
Natl  Acad  Sci  USA  2010;  107:  2669-2674. 


12  Hammamieh  R,  Chakraborty  N,  De  Lima  TC,  Meyerhoff  J,  Gautam  A,  Muhie  S  et  al. 
Murine  model  of  repeated  exposures  to  conspecific  trained  aggressors 
simulates  features  of  post-traumatic  stress  disorder.  Behav  Brain  Res  2012;  235: 
55-66. 

13  Chakraborty  N,  Meyerhoff  J,  Gautam  A,  Muhie  S,  Jibitu  M,  De  Lima  TC  et  al.  Gene 
and  stress  history  interplay  in  emergence  of  PTSD-like  features.  Behav  Brain  Res 
2015;  292:  266-277. 

14  Muhie  S,  Gautam  A,  Meyerhoff  J,  Chakraborty  N,  Hammamieh  R,  Jett  M.  Brain 
transcriptome  profiles  in  mouse  model  simulating  features  of  post-traumatic 
stress  disorder.  Mol  Brain  2015;  8:  14. 

15  Magni  P,  Ferrazzi  F,  Sacchi  L,  Bellazzi  R.  TimeClust:  a  clustering  tool  for  gene 
expression  time  series.  Bioinformatics  2008;  24:  430-432. 

16  Xia  J,  Gill  EE,  Hancock  RE.  NetworkAnalyst  for  statistical,  visual  and  network-based 
meta-analysis  of  gene  expression  data.  Nat  Protoc  2015;  10:  823-844. 

17  Kramer  A,  Green  J,  Pollard  J  Jr,  Tugendreich  S.  Causal  analysis  approaches  in 
Ingenuity  Pathway  Analysis.  Bioinformatics  2014;  30:  523-530. 

18  Liu  Q,  Tan  Y,  Huang  T,  Ding  G,  Tu  Z,  Liu  L  et  al.  TF-centered  downstream  gene  set 
enrichment  analysis:  Inference  of  causal  regulators  by  integrating  TF-DNA  inter¬ 
actions  and  protein  post-translational  modifications  information.  BMC  Bioinfor¬ 
matics  2010;  11(Suppl  11):  S5. 

19  Lahti  L,  Haugas  M,  Tikker  L,  Airavaara  M,  Voutilainen  MH,  Anttila  J  et  al.  Differ¬ 
entiation  and  molecular  heterogeneity  of  inhibitory  and  excitatory  neurons 
associated  with  midbrain  dopaminergic  nuclei.  Development  2016;  143:  516-529. 

20  Kohman  RA,  Rhodes  JS.  Neurogenesis,  inflammation  and  behavior.  Brain  Behav 
Immun  2013;  27:  22-32. 

21  Jessberger  S,  Clark  RE,  Broadbent  NJ,  Clemenson  GD,  Consiglio  A,  Lie  DC  et  al. 
Dentate  gyrus-specific  knockdown  of  adult  neurogenesis  impairs  spatial  and 
object  recognition  memory  in  adult  rats.  Learn  Memory  2009;  16:  147-154. 

22  Bremner  JD,  Elzinga  B,  Schmahl  C,  Vermetten  E.  Structural  and  functional  plasticity 
of  the  human  brain  in  posttraumatic  stress  disorder.  Prog  Brain  Res  2008;  167: 
171-186. 

23  Miller  MM,  McEwen  BS.  Establishing  an  agenda  for  translational  research  on  PTSD. 
Ann  NY  Acad  Sci  2006;  1071:  294-312. 

24  Shin  LM,  Rauch  SL,  Pitman  RK.  Amygdala,  medial  prefrontal  cortex,  and  hippo¬ 
campal  function  in  PTSD.  Ann  NY  Acad  Sci  2006;  1071:  67-79. 

25  Michopoulos  V,  Powers  A,  Gillespie  CF,  Ressler  KJ,  Jovanovic  T.  Inflammation  in 
fear-  and  anxiety-based  disorders:  PTSD,  GAD,  and  beyond.  Neuropsycho¬ 
pharmacology  201 7;  42:  254-270. 

26  Kirby  ED,  Friedman  AR,  Covarrubias  D,  Ying  C,  Sun  WG,  Goosens  KA  et  al.  Baso- 
lateral  amygdala  regulation  of  adult  hippocampal  neurogenesis  and  fear-related 
activation  of  newborn  neurons.  Mol  Psychiatry  2012;  17:  527-536. 

27  Van  der  Borght  K,  Havekes  R,  Bos  T,  Eggen  BJ,  Van  der  Zee  EA.  Exercise  improves 
memory  acquisition  and  retrieval  in  the  Y-maze  task:  relationship  with  hippo¬ 
campal  neurogenesis.  Behav  Neurosci  2007;  121:  324-334. 

28  Adeosun  SO,  Hou  X,  Zheng  B,  Stockmeier  C,  Ou  X,  Paul  I  et  al.  Cognitive  deficits 
and  disruption  of  neurogenesis  in  a  mouse  model  of  apolipoprotein  E4  domain 
interaction.  J  Biol  Chem  2014;  289:  2946-2959. 

29  Graciarena  M,  Depino  AM,  Pitossi  FJ.  Prenatal  inflammation  impairs  adult  neu¬ 
rogenesis  and  memory  related  behavior  through  persistent  hippocampal  TGF 
beta(1)  downregulation.  Brain  Behav  Immun  2010;  24:  1301-1309. 

30  Sakalem  ME,  Seidenbecher  T,  Zhang  M,  Saffari  R,  Kravchenko  M,  Wordemann  S 
et  al.  Environmental  enrichment  and  physical  exercise  revert  behavioral  and 
electrophysiological  impairments  caused  by  reduced  adult  neurogenesis.  Hippo¬ 
campus  2017;  27:  36-51. 

31  Gautam  A,  D'Arpa  P,  Donohue  DE,  Muhie  S,  Chakraborty  N,  Luke  BT  et  al.  Acute 
and  chronic  plasma  metabolomic  and  liver  transcriptomic  stress  effects  in  a 
mouse  model  with  features  of  post-traumatic  stress  disorder.  PLoS  ONE  2015;  10: 
eOl  17092. 

32  Cho  JH,  Lee  I,  Hammamieh  R,  Wang  K,  Baxter  D,  Scherler  K  et  al.  Molecular 
evidence  of  stress-induced  acute  heart  injury  in  a  mouse  model  simulating 
posttraumatic  stress  disorder.  Proc  Natl  Acad  Sci  USA  2014;  111:  3188-3193. 

33  Widera  D,  Mikenberg  J,  Kaltschmidt  B,  Kaltschmidt  C.  Potential  role  of  NF-kappa  B 
in  adult  neural  stem  cells:  the  underrated  steersman?  Int  J  Dev  Neurosci  2006;  24: 
91-102. 

34  Baron  R,  Nemirovsky  A,  Harpaz  I,  Cohen  H,  Owens  T,  Monsonego  A.  IFN-gamma 
enhances  neurogenesis  in  wild-type  mice  and  in  a  mouse  model  of  Alzheimer's 
disease.  FASEB  J  2008;  22:  2843-2852. 

35  Erta  M,  Quintana  A,  Hidalgo  J.  Interleukin-6,  a  major  cytokine  in  the  central 
nervous  system.  Int  J  Biol  Sci  2012;  8:  1254-1266. 

36  Bowen  KK,  Dempsey  RJ,  Vemuganti  R.  Adult  interleukin-6  knockout  mice  show 
compromised  neurogenesis.  Neuroreport  2011;  22:  126-130. 

37  Baker  DG,  Nievergelt  CM,  O'Connor  DT.  Biomarkers  of  PTSD:  neuropeptides  and 
immune  signaling.  Neuropharmacology  2012;  62:  663-673. 

38  Yirmiya  R,  Goshen  I.  Immune  modulation  of  learning,  memory,  neural  plasticity 
and  neurogenesis.  Brain  Behav  Immun  2011;  25:  181-213. 


Translational  Psychiatry  (2017),  1  -9 


Stress-induced  neuroinflammation 
S  Muhie  et  al 


39  Wood  SK,  Wood  CS,  Lombard  CM,  Lee  CS,  Zhang  XY,  Finnell  JE  et  al.  Inflammatory 
factors  mediate  vulnerability  to  a  social  stress-induced  depressive-like  phenotype 
in  passive  coping  rats.  Biol  Psychiatry  2015;  78:  38-48. 

40  Groer  MW,  Kane  B,  Williams  SN,  Duffy  A.  Relationship  of  PTSD  symptoms  with 
combat  exposure,  stress,  and  inflammation  in  American  soldiers.  Biol  Res  Nurs 
2015;  17:  303-310. 

41  Wu  MD,  Hein  AM,  Moravan  MJ,  Shaftel  SS,  Olschowka  JA,  O'Banion  MK.  Adult 
murine  hippocampal  neurogenesis  is  inhibited  by  sustained  IL-1  beta  and  not 
rescued  by  voluntary  running.  Brain  Behav  Immun  2012;  26:  292-300. 

42  Belarbi  K,  Rosi  S.  Modulation  of  adult-born  neurons  in  the  inflamed  hippocampus. 
Front  Cell  Neurosci  2013;  7:  145. 

43  Parihar  VK,  Hattiangady  B,  Shuai  B,  Shetty  AK.  Mood  and  memory  deficits  in  a 
model  of  Gulf  War  illness  are  linked  with  reduced  neurogenesis,  partial  neuron 
loss,  and  mild  inflammation  in  the  hippocampus.  Neuropsychopharmacology 
2013;  38:  2348-2362. 

44  Zonis  S,  Pechnick  RN,  Ljubimov  VA,  Mahgerefteh  M,  Wawrowsky  K,  Michelsen  KS 
et  al.  Chronic  intestinal  inflammation  alters  hippocampal  neurogenesis.  J  Neu¬ 
roinflammation  2015;  12:  65. 

45  Ho  N,  Sommers  MS,  Lucki  I.  Effects  of  diabetes  on  hippocampal  neurogenesis: 
links  to  cognition  and  depression.  Neurosci  Biobehav  Rev  2013;  37:  1346-1362. 

46  Passos  1C,  Vasconcelos-Moreno  MP,  Costa  LG,  Kunz  M,  Brietzke  E,  Quevedo  J  et  al. 
Inflammatory  markers  in  post-traumatic  stress  disorder:  a  systematic  review, 
meta-analysis,  and  meta-regression.  Lancet  Psychiatry  2015;  2:  1002-1012. 

47  Gill  JM,  Saligan  L,  Woods  S,  Page  G.  PTSD  is  associated  with  an  excess  of 
inflammatory  immune  activities.  Perspect  Psychiatr  Care  2009;  45:  262-277. 

48  Pace  TW,  Wingenfeld  K,  Schmidt  I,  Meinlschmidt  G,  Hellhammer  DH, 
Heim  CM.  Increased  peripheral  NF-kappaB  pathway  activity  in  women  with 
childhood  abuse-related  posttraumatic  stress  disorder.  Brain  Behav  Immun  2012; 
26:  13-17. 

49  Lindqvist  D,  Wolkowitz  OM,  Mellon  S,  Yehuda  R,  Flory  JD,  Henn-Haase  C  et  al. 
Proinflammatory  milieu  in  combat-related  PTSD  is  independent  of  depression 
and  early  life  stress.  Brain  Behav  Immun  2014;  42:  81-88. 

50  Rohleder  N,  Karl  A.  Role  of  endocrine  and  inflammatory  alterations  in  comorbid 
somatic  diseases  of  post-traumatic  stress  disorder.  Minerva  Endocrinol  2006;  31: 
273-288. 

51  Voiko  CD,  Regidor  PA,  Rohr  UD.  Model  approach  for  stress  induced  steroidal 
hormone  cascade  changes  in  severe  mental  diseases.  Horm  Mol  Biol  Clin  Investig 
2016;  25:  157-170. 

52  Sharp  TJ.  The  prevalence  of  post-traumatic  stress  disorder  in  chronic  pain 
patients.  Curr  Pain  Headache  Rep  2004;  8:  111-115. 


53  Beckham  JC,  Crawford  AL,  Feldman  ME,  Kirby  AC,  Hertzberg  MA,  Davidson  JR  et  al. 
Chronic  posttraumatic  stress  disorder  and  chronic  pain  in  Vietnam  combat 
veterans.  J  Psychosom  Res  1997;  43:  379-389. 

54  Sharp  TJ,  Harvey  AG.  Chronic  pain  and  posttraumatic  stress  disorder:  mutual 
maintenance?  Clin  Psychol  Rev  2001;  21:  857-877. 

55  Brennstuhl  MJ,  Tarquinio  C,  Montel  S.  Chronic  pain  and  PTSD:  evolving  views  on 
their  comorbidity.  Perspect  Psychiatr  Care  2015;  51:  295-304. 

56  Asmundson  GJ,  Coons  MJ,  Taylor  S,  Katz  J.  PTSD  and  the  experience  of  pain: 
research  and  clinical  implications  of  shared  vulnerability  and  mutual 
maintenance  models.  Can  J  Psychiatry  2002;  47:  930-937. 

57  Plantinga  L,  Bremner  JD,  Miller  AH,  Jones  DP,  Veledar  E,  Goldberg  J  et  al.  Asso¬ 
ciation  between  posttraumatic  stress  disorder  and  inflammation:  a  twin  study. 
Brain  Behav  Immun  2013;  30:  125-132. 

58  Wilson  CB,  McLaughlin  LD,  Nair  A,  Ebenezer  PJ,  Dange  R,  Francis  J.  Inflammation 
and  oxidative  stress  are  elevated  in  the  brain,  blood,  and  adrenal  glands  during 
the  progression  of  post-traumatic  stress  disorder  in  a  predator  exposure 
animal  model.  PLoS  ONE  2013;  8:  e76146. 

59  Brudey  C,  Park  J,  Wiaderkiewicz  J,  Kobayashi  I,  Mellman  TA,  Marvar  PJ.  Autonomic 
and  inflammatory  consequences  of  posttraumatic  stress  disorder  and  the  link  to 
cardiovascular  disease.  Am  J  Physiol  Regul  Integr  Comp  Physiol  2015;  309: 
R315-R321. 

60  O'Donovan  A,  Chao  LL,  Paulson  J,  Samuelson  KW,  Shigenaga  JK,  Grunfeld  C  et  al. 
Altered  inflammatory  activity  associated  with  reduced  hippocampal  volume  and 
more  severe  posttraumatic  stress  symptoms  in  Gulf  War  veterans.  Psychoneur¬ 
oendocrinology  2015;  51:  557-566. 

61  Yang  R,  Daigle  BJ  Jr.,  Muhie  SY,  Hammamieh  R,  Jett  M,  Petzold  L  et  al.  Core 
modular  blood  and  brain  biomarkers  in  social  defeat  mouse  model  for  post 
traumatic  stress  disorder.  BMC  Syst  Biol  2013;  7:  80. 


This  work  is  licensed  under  a  Creative  Commons  Attribution- 
NonCommercial-ShareAlike  4.0  International  License.  The  images  or 
other  third  party  material  in  this  article  are  included  in  the  article’s  Creative  Commons 
license,  unless  indicated  otherwise  in  the  credit  line;  if  the  material  is  not  included  under 
the  Creative  Commons  license,  users  will  need  to  obtain  permission  from  the  license 
holder  to  reproduce  the  material.  To  view  a  copy  of  this  license,  visit  http:// 
creativecommons.org/licenses/by-nc-sa/4.0/ 


©  The  Author(s)  2017 


Supplementary  Information  accompanies  the  paper  on  the  Translational  Psychiatry  website  (http://www.nature.com/tp) 


Translational  Psychiatry  (2017),  1  -9 


